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Abstract 

We construct a fully back-reacted holographic dual of a four-dimensional field theory 
which exhibits chiral symmetry breaking. Two possible models are considered by studying 
the effects of a five-dimensional field, dual to the qq operator. One model has smooth 
geometry at all radii and the other dynamically generates a cutoff at finite radius. Both 
of these models satisfy Einstein's field equations. The second model has only three free 
parameters, as in QCD, and we show that this gives phenomenologically consistent results. 
We also discuss the possibility that in order to obtain linear confinement from a back- 
reacted model it may be necessary to consider the condensate of a dimension two operator. 
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1 Introduction 



Strong interactions are described in the standard model (SM) by an SU (3) gauge theory known 
as quantum chromodynamics (QCD) [T]. Quarks, which are basic constituents in the SM, are in 
the fundamental representation of the SU (3) group. They interact with each other by coupling 
to the gluons, the vector gauge bosons transforming as the adjoint representation of SU{3). 
Since the gauge group is non-Abelian the gluons have direct self-interactions (this feature is 
closely related to the origin of life). It is these self- interactions that cause a negative beta 
function, for the running coupling constant as{fi) [21 [3]. This is a unique feature among 

four-dimensional renormalizable gauge theories. Because < 0, the coupling constant as{fi) 
decreases at short distances (UV region) and it is this anti-screening effect, known as asymptotic 
freedom, which makes the UV region perturbatively understood. 

At low energies (IR region), solving the theory becomes extremely difficult. As the coupling 
constant as(/i) grows in the IR, perturbative methods are no longer applicable. We are currently 
unable to solve from first principle the low energy dynamics of QCD. One can only study the low 
energy features of strong interactions from experimental data and construct effective quantum 
field theories to characterize the low energy features of QCD, such as dynamically generated 
spontaneous symmetry breaking [1] . It was recently shown in ref. [5] that the gap equation can 
be derived as the minimal condition of the effective Higgs potential generated dynamically and 
the lightest scalar mesons as the parters of pseudoscalar mesons play the role of the composite 
Higgs bosons to realize the dynamical spontaneous symmetry breaking. As a consequence, the 
mass spectrum of scalar and pseudoscalar mesons and their mixings as well as the relevant low 
energy parameters were predicted |5] to be consistent with the experimental data by using the 
same number of input parameters as in QCD. At low energies, quarks are bound to form hadrons, 
which are color singlets composed of quarks and gluons. Isolated quarks and gluons appear not 
to exist. This so called color confinement phenomenon has been a well known conjecture since 
the development of QCD and has not yet been proved. 

The details of hadronization are still unclear. In the massless limit, the theory has only one 
parameter, the coupling constant. If QCD is the fundamental theory of the strong interactions, 
it must have an intrinsic scale to account for the short-range feature of the strong interactions, 
even if all the fields are massless. This is a pure quantum effect and has the consequence that the 
lightest hadron is massive. In the QCD Lagrangian, one can see that quark masses act as sources 
of the operators q{x)q{x). This means that if operators q{x)q{x) are allowed, quarks must be 
massive. In other words, q{x)q{x) must break some symmetry in the massless limit. It is not 
difficult to see that this symmetry in the massless limit is the global SU (Nj) i SU {Nf)ji chiral 
symmetry. Therefore, it is natural to expect that the formation of the condensate {q{x)q{x)) 
will dynamically break the chiral symmetry. None of the non-perturbative features mentioned 
above has yet been justified from first principles. 

In p] 't Hooft proposed the idea that one might be able to solve QCD by first studying an 
SU{N) gauge theory in the limit 0. One can then perform a expansion and study 

the theory perturbatively. In this limit, one can argue that properties such as color confinement 
and dynamical chiral symmetry breaking still exist. Furthermore, the spectra are given by an 
infinite number of meson and glueball resonances. Mesons are stable and the decay amplitudes 
are 0(1/ a/ZV). The large theory can be treated as an effective theory of hadrons and glueballs. 
Since the intrinsic scale Aqcd is independent of when 1/A^ — > 0, the mass of mesons does 
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not change with A^. It can be shown that current-current correlation functions can be written 
as an infinite sum over meson resonances. The poles and the corresponding residues give the 
mass and decay constant of mesons, if one knows how to calculate the two-point function. 

The duality between gravity and gauge theories conjectured by Maldacena [7] and further 
developed in [SI [1] has shed new light on the problem of strongly coupled gauge theories. 
Although the original conjecture in [7] is the mathematical correspondence between the low 
energy approximation of type IIB string theory on AdS^ and N" = 4 U (N) SYM for 

large N in four dimensions, the holographic principle can in theory be applied to general gauge 
theories and their gravity duals. The idea goes as follows. Consider a manifold Ai whose 
limit in infinity is X with X being an Einstein manifold and W being compact. After 

compactification on X, one can build a manifold ^ with boundary S such that the metric on 
X has a double pole near the boundary S. Then the isometry group of X acts as the conformal 
group of S. The duality is the holographic correspondence that for a conformal field theory 
(CFT) on S, one can find a gravity theory on X such that the generating functional of the 
connected Greens functions in CFT on S is equivalent to the action of the gravity theory on A?, 
whose fields on S are identified to be the sources of the operators in CFT (see [TU] for recent 
work exploring generalisations of this relationship). The mass of a p-form field on A" is 
quantized and related to the conformal dimension A of the operator on S it corresponds to. 
The relation is [9] = {A — p){A + p — d) where d is the dimension of S. The anomalous 

dimension of the operator on S corresponds to the quantum corrections on m<^. That is, the 
quantum corrections on m<^, (5m<^, should satisfy 

r /A d^ '^ddtn 

bm^ = (Ao 
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where is the anomalous dimension of O^^^ and we use the subscript to denote the un- 
renormalized quantities. If the anomalous dimension of an operator vanishes, the operator must 
be a generator of some symmetry algebra on S. In this case, the must be protected from 
corrections by some symmetry on X. For example, the dimension 3 operators q'j'^t"'q are con- 
served currents in QCD. Their holographic dual are 1-form fields whose masses vanish from 
the relation mentioned above. must be vector gauge fields on X and remain massless by 
gauge symmetry. 

QCD is clearly not a conformal field theory. However, due to the asymptotic freedom, it 
approaches to the conformal limit in UV regions. At high energies, one may neglect quark 
masses and the chiral symmetry is restored. In these regions, when approaching the fixed point, 
which is suggested to exist [TTl [T2] . QCD becomes a conformal field theory and the holographic 
recipe is applicable. Many efforts have been made in constructing a holographic dual of QCD. 
See p!3l ] and references therein for various approaches. Usually in these models, the background 
manifold is chosen to be AdS, which means one did not consider the effects of the condensates. 
This, however, is not logically consistent. In QCD, the relevant and marginal operators that can 
form condensates without breaking Lorentz and gauge symmetries at low energies are F^^F"'^'^ 
(dimension 4) and qq (dimension 3). The dimension 2 operator condensate (A^A°^)mm[27] 
will be discussed later. These condensates can be used to parameterize the non-perturbative 
observables and their effects should not be neglected. The holographic dual of these operators 
are scalar fields, acting as sources for these operators on the boundary of X. Including the 
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effects of these condensates means we are introducing these scalar fields in the bulk, which in 
turn will deform the AdS manifold. In this article, we consider the effects of these scalar fields 
and solve the graviton-tachyon system to give a correct deformed AdS metric and perform a 
consistent analysis. Adopting this bottom-up approach, one is running a risk that the model 
built might not have anything to do with the underlying theory. Nevertheless, if the model can 
describe the distinctive features of QCD and accurately predict the non-perturbative quantities 
such as mass spectrum and decay constants of hadrons, then it deserves one's attention and 
might be able to give us some clues about the structure of the underlying theory if one believes 
that QCD has a holographic dual theory. 



2 Effects of The Condensates 

Before constructing a phenomenological model, we will first probe the general effects of a bulk 
field (/)(a;^, y) on the metric, at the same time setting up the notation. We start with the following 
five- dimensional effective action 

S = l_ d'x^i-R + 1(90)2 + Vi<P)) , (1) 

where R is the five- dimensional Ricci scalar. 

With time-reversal and parity symmetries, one can show that the static metric which re- 
spects Poincar symmetry on four- dimensional coordinates x'^, the one we are interested, has the 
following form 

ds^ = e-^^^y^x^dx^" - dy^ . (2) 

As we stressed before, we are building a holographic model based on our knowledge of QCD, the 
potential V{(j)) is unknown and will be determined once the specific form of y) is chosen to 
describe some QCD condensate. Without matter fields, X is an Einstein manifold. Demanding 
X to be AdS^ in that case, it is easy to show that the warp factor A{y) = y and the constant 
potential V equals A = 12 in our units. (With arbitrary dimension D of X, A = {D — 1){D — 2).) 
The equations of motion for Eqn. ([T]) are 

IdPQi-R + ^^M^a^V + V{<P)] + RpQ - ^dpcj>dQ<P = , (3) 



and 



which gives 



6A"(y)-12A'2(y) + y(0)-i0'2 = O, (5) 
12A'\y) - Vi<P) = ^<P'' , (6) 
0"(y)-4A'(,)0'(,) + ^ = O. (7) 
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Here we consider the particular solution such that (j) is independent of x^, since we demand 
QCD condensates are translationally invariant. The equations can be solved using the method 
in [28]. Eqn.([5]) and Eqn.([6]) gives 

qX' - 0'' = , (8) 

and 

3A" - 12A'2 + V{(t>) = . (9) 

Note that care must be taken here as the factor of 6 in Eqn.([H]) becomes a 3 in the case of 
a complex field which is used in section [31 This factor also affects the following equations in 
this section. As we know that without the scalar field 0, = A and A {y) = 1 (maximally 
symmetric). Therefore, it is natural to expect that A (y) = W{(j){y)) after the effect of (j) is 
considered. Eqn.(IH]) tells us that 6^ = (f)' since A" = ^^4> ■ Eqn.(l9]) can be taken as the 
definition of ^(0), that it, 1/(0) = UW"^ - 18(^)2. It is trivial to check that Eqn.© will 
automatically be satisfied: 

In the UV region {y — oo), the effects of condensates are less and less important so we should 
expect (piy) — as ?/ — — oo. In fact, from Eqn.Q one can see that when y — » — oo, keeping 
only the first nontrivial term of V{(j)) (i.e., up to 0(0^)), the linearized equation of motion for 
(f){y) gives (j){y) ~ e"^. For 0(?/) = ae"^, = nip. And it is easy to show that W{(j)) = y|0^ + 1 
so that A{y) = |^e^"^ + y. Therefore, we arrive at an effective action whose potential V{(j)) 
equals ^(0) = 12 + (2n — ^)0^ + 0^- This can be easily generalized to the case with n scalar 
fields (pi. For 0j = aic"-^^, the potential V^(0i) is 

Vi^i) = 12 + 5^(2n, - ^nM + ^ E ^^M^j + J^H 

i<j 



and the warp factor A{y) is 



A{y)=y + ^Y.^y^'"- 



The constant term guarantees that without 0j, X goes back to AdS^. The mass terms only 
depend on Ui. The model we build in this way should be the holographic dual of a four- 
dimensional theory. Note that this example for is an oversimplification and a true dual field 
will include at least two terms. Each field 0, on the UV boundary should act as the source of 
some dimension Aj operator on the four- dimensional theory. According to m^, = A, (A, — 4). 
Therefore, we have rij = Aj or 4 — Aj and the meaning of rii is clear now. Also from the mass 
terms one can see that the condensates formed from relevant operators correspond to tachyonic 
scalar fields on the holographic theory. As long as one demands (i) the warp factor A{y) = y (i.e., 
AdS^) when matter fields (f)i{y) = and (ii) (j)i{y) ~ aie^'^ (or e*^^"^^)^) in the UV {y — > — oo) 
in order to do the correspondence, the first two terms of the potential V{(f)i) are determined, 
independent of the specific forms of 0j. From this example one might think that the coupling 
constants of the self- interaction terms only depend on n^. This, however, is in general not true. 

In the following section, we will use the method discussed in this section to construct a 
phenomenological model for mesons. 
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3 The Phenomenological Model 



We wish now to consider the phenomenological implications of including a back-reaction of a 
five-dimensional field on the originally AdS geometry. We must first decide which fields are 
relevant to the theory in question. In QCD there are four relevant operators which we may 
consider given by qlq^ (lL,Rl^t°'(lL,R, F"^ and FF. These have been discussed in the hterature 
in [131 [H]- In this paper we concentrate on the complex scalar, dimension three operator 
to study its effect on the geometry and subsequently the phenomenology of the model. The 
inclusion of the source for a dimension four operator in addition is no more complicated but it 
does introduce another free parameter, the condensate of F"^. Though this operator is clearly 
important, we can argue that its effect will be of a similar order of magnitude to that of the 
dimension three operator. 

We define the source for the qlq^ operator as a matrix valued field X{x^, y) proportional to 
the unit two by two matrix. This can be extended to three flavours but for now we discuss only 
two. The action is given by: 



S = fd^x^Tr 



.R+\DX\^ - V{X) - ^^{Fl + Fl) 



(10) 



V{X) = -?>\x\^ + vux) (11) 

where the self-interaction terms are to be determined by a consistent back-reaction and the 
vector and axial vector gauge fields are sources of the vector and axial vector currents in the 
four-dimensional theory. 

We must determine a consistent form for X which will provide a sound phenomenological 
model. There are several distinct routes one can take at this point, each of which has different 
repercussions for the model. The first option is to make an ansatz for X which is polynomial in 
z{= e^) with a finite number of terms. The second option is to construct a function of z which 
has an infinite series expansion but has no poles. The third option is to define a function which 
is singular at some specific point z in the geometry. Each of these will be discussed in turn. 

We know that the AdS behavior of the field X, corresponding to a dimension three scalar 
operator must go like: 

X=(^z + |.^)W (12) 

where the coefficients can be shown to correspond to a quark mass and bilinear condensate 
respectively. If we wish to construct a polynomial expression with a finite number of terms, 
greater than two, we will have to introduce more free parameters to our model which will 
necessarily be fixed by experimental data. These additional terms will make the model less 
predictive, so it appears that if we want a finite polynomial expression, it is most natural to 
keep the AdS form of the field X. Indeed we will show that this can be done consistently while 
finding the back-reaction on the geometry. 

The second option may have significant implications in the context of linear confinement 
where a specific IR behavior for the metric and dilaton have been discussed in [22]. Finding a 
back-reacted solution which gives this is an important goal and is discussed more in section HI 

The third option provides an interesting prospect, by which the space is truncated smoothly 
without the need of an IR boundary put in by hand, though the interpretation of the additional 
field may itself come from some distribution of branes in the IR. The benefit of such a choice 
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Figure 1: The potential in the five-dimensional action for both the non-back- reacted (dashed 
line) and back-reacted (full line) theories. 



may be that the number of free parameters of the model is reduced, by making the value of the 
IR cutoff depend explicitly on the value of the quark bilinear condensate, as is expected from 
QCD. The question however is, what function should one choose for the X field. We show in 
section 13.21 that a simple guess gives a consistent set of predictions. 

3.1 Model I 

From Eqn(l8]) we can calculate the effect of the tachyonic field given by Eqn( |T2l) on the metric 
and find that the back-reacted function A is given by: 

A{z) = log{z) + 4 ^ 6 (13) 

V ; ^ 24 16 24 ^ ^ 

Figure [1] shows the potential in this back-reacted theory compared to the potential without the 
back-reacton. The potential goes roughly as for large X. The inclusion of the back-reaction 
now means that the mass of the vector mesons will be a function of the quark mass and quark 
bilinear condensate. The question is, how much of an effect this will be. One expects in QCD 
that the mass of the vector mesons should be determined from the value of Aqcd and should 
not be affected significantly by the quark parameters. 

Just as in [13] we choose to fix the four free parameters in this theory using Nc, fn, and 
Mp. A simple check shows that, as expected, the addition of the back-reaction does not affect the 
proof that the Gell-Mann-Oakes-Renner relation holds in this model. This immediately allows 
us to fix the quarks bilinear condensate as a function of the quark mass, f-j^ and (errors 
will be of the order and therefore negligible in any reasonable model). The five-dimensional 
coupling is fixed from the vector current two-point function to be 

9l = ^^ (14) 

and for the phenomenological model we pick = 3. This leaves us with two free parameters, 
niq and zm, the IR cutoff. 
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Figure 2: ip{e) as a function of zjr and nig, showing the contour which corresponds to the 
physical boundary value for this field. 



To fix zjFt we use the position of the pole in the vector current correlator, provided by the 
eigenvalue of its dual five-dimensional field. The equation of motion for an arbitrary component 
of the vector field, ip{z), is given by Eqn ffTSl) . 

12e^^(^'''''^^^'^''''^^<''hiP{z)ml- {Sa^z^^ + Sam^z"" + my + 12) ij'{z) + 12z^p"{z) = , (15) 



where we have the constraint that 



a 



2mn 



(16) 



The boundary conditions for the field if) are given by a vanishing derivative at the IR-brane and 
a vanishing value in the UV (at 2 — > 0). In order to calculate the p-meson decay constant, we 
must be careful of the normalization, however, the normalization will not affect the position of 
the zeros of the function in the UV and we are therefore free to choose ipizm) = 1. Having 
fixed the IR conditions of ip completely, we can vary zjpi and rrig to find which values give the 
correct z behavior. By plotting the value of ip{e) as a function of the two free parameters 
on a scale such that all we can see is the switch from positive to negative it is clear that 
this is almost independent of the quark mass, for a reasonable range of values. This result is 
as expected and tells us that the vector spectrum is determined largely through Aqcd and not 
through rriq. 

In order to fix the two parameters completely we now turn to the pion decay constant, 
calculated from the axial vector equation of motion with the pole in the propagator set to zero. 
Note that the axial vector equation of motion depends on the quark mass and condensate even 
neglecting the back reaction. The pion decay constant is given by: 

1 d,Aiq = 0,z) 



f 



9l 



(17) 
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Where the solution of the axial field A is the solution of 

4:871^ zA{z) [az^ + m^f + (Sa^^^ + ^am^z^ + m^^^ ^ 12) A'{z) - 12zA"{z) = . (18) 

This time we plot the value of fn,mod— fn,exp as a function of the two free parameters, where fn,mod 
corresponds to the value calculated in this model and fn,exp corresponds to the experimental 
value. By combining this plot with the plot in figure [2] we can fix the values of the two free 
parameters precisely. These are found to be rriq = 2AMeV and cr = (323.5Mel^)^. We fix the 
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Figure 3: Intersection of the physical contours for both and ip^e) as a function of the two 
free parameters zjji and rriq. The point of intersection fixes the values to be used to calculate 
the observable predictions. 

remaining free parameters of the theory from the intercept values of the two graphs and are 
now ready to calculate the remaining observable quantities, precisely as in [13]. The results of 
this analysis are given in table [U 

We see that the effect of the back-reaction in this case is very small. We have shown that 
we can construct a phenomenologically consistent model, including back-reaction and obtain 
results in very close agreement with those without the back-reaction. The reason for this is 
simple. The effect of the back-reaction on the geometry will only become significant for energies 
of the order of era. With the inclusion of the IR-brane, the geometry is cutoff just below the 
region where these effects would become significant and so the effects of the quark dynamics 
are shown to be negligible. 

3.2 Model II 

As stated in the previous section, it would be appealing to remove zjr as a free parameter and 
instead to generate it dynamically via a condensate. We examine one possible function with 
which we can do this and show that the results are promising. 
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Table 1: Results of the non-back-reacted, non-singular back-reacted (Model I) and singular 
back-reacted (Model II) geometries. The percentage error is taken from the central value of the 
experimental measurement. * indicates that this observable was used to fit the data. 



Motivated by supergravity solutions we make a simple choice for X, of the form 



where the limiting behavior is for small z. The space is therefore cutoff at zir = 

0-3 

Again we use the Gell-Mann-Oakes-Renner relation to constrain the relationship between 
rrig and a leaving us with one free parameter which we can fix using the pion decay constant. 
The values are found to be = 2.53 MeV and a = (320.3Mel^)^. 

We perform the same analysis using this model as we did in section 13.11 The results of 
this analysis are given in table [H This singular back-reacted geometry gives good results and 
the important point is that the singular behavior does not destroy the original phenomenology. 
There are clearly many other sensible choices for X which are both singular and give the correct 
UV behavior, however we show that with a sensible choice, such a model is phenomenologically 
consistent. 

We can study the potential in the case of this singular back-reacted geometry and compare 
it to the original potential in the non-back-reacted case. This is shown in figure HI We see that 
while the original potential is unbounded, the potential which classically gives us the singular 
field, X, gives a bounded potential with a minimum. The position of the minimum is a function 
of rriq and a but its existence is not dependent on these values. Though we treat the five- 
dimensional theory classically, from a quantum mechanical perspective, the potential in this 
singular, back-reacted geometry would have important consequences. 



4 Discussion and Conclusion 

In this paper we consider the back-reaction of the scalar on the metric in the holographic QCD 
model. In the hard wall model, the meson spectrum is generated by the IR brane. For pure 
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Figure 4: The potential in the five-dimensional action for both the non-back-reacted (dashed 
line) and back- reacted (full line), Model II, theories. 



AdS without the hard wall, the five-dimensional model is scale invariant and all the mesons 
are massless. Breaking the scale invariance is an essential mechanism in order to generate the 
meson spectrum. The mass gap Aqcd is dual to the IR cut-off in this model. After considering 
the back-reaction, we find that the model still successfully describes the spectrum and decay 
constants of ground state mesons. 

As for the resonances, in [22] it was shown that one can get linear confinement by including 
a quadratic dilaton $ ~ 2;^ in the pure AdS background (see also [26] for a phenomenological 
treatment including a UV cutoff). The peculiar dependence of the dilaton is far from clear. 
If one considers the deformed AdS metric without introducing a dilaton: 

ds^ = l^fdx.dx" -dz^), 
z^ 

it can be shown that with /(z) ~ (cosh(a2r^) -|- 6sinh(a2^))'^ and a 7^ 0, one breaks scale 
invariance and gets linear confinement, even without the presence of the IR brane. However, 
the deformed AdS metric given by the back-reaction of the quark condensate is not compatible 
with this. Linear confinement can not be realized in this phenomenological model. From the 
non-polynomial form of ]{^z\ the first non-trivial term is quadratic in . This suggests that a 
non-local dimension two condensate should be formed. In QCD, it has been investigated [27| 
that the non-zero value for the minimum of the vector potential condensate A? ^ which is defined 
to be J (i^x/l^(x) is possible. This dimension two condensate can be shown to be gauge invariant 
but non-local. Therefore, its effect should be relevant to linear confinement. 
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